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COMPARISON OF PREDICTION METHODS AND STUDIES OF 
RELAXATION IN HYPERSONIC TURBULENT 


NOZZLE -WALL BOUNDARY LAYERS* 

By Dennis M. Bushnell, Charles B. Johnson, 

William D. Harvey, and William V. Feller 
Langley Research Center 

SUMMARY 

Recent turbulent boundary-layer measurements on axisymmetric nozzle walls at 
Mach numbers of 6, 8, and 19 are presented and compared with the predictions of two 
theoretical methods. One of the methods was of the integral type; the other was a finite- 
difference solution to the governing differential equations. The integral method gave 
good predictions of the data at Mach 6 and 8 and a reasonable prediction at Mach 19. 
Results for the finite -difference approach were available only at Mach 8, for which good 
agreement with data was obtained. 

A survey of boundary -layer data indicates that the variation of total temperature 
with velocity for flat-plate type of flows can be represented within the scatter of the data 
by the linear Crocco relation, whereas for nozzle-wall flows, this relation is generally 
more nearly quadratic. A finite -difference solution for the Mach 8 nozzle-wall boundary 
layer resulted in a quadratic -type temperature-velocity relation in the strong favorable 
pressure -gradient region downstream of the throat. However, near the nozzle exit the 
theoretical solution tended to approach the linear variation of Crocco. Apparently, the 
observed quadratic variation in nozzle boundary layers can be at least partly accounted 
for by the pressure-gradient history of the flow. 

Profile measurements along a straight pipe extension downstream of the nozzle 
exit of the Mach 6 tunnel indicate that at a distance of the order of 60 boundary-layer 
thicknesses, the flow tended to relax toward the linear temperature-velocity variation. 

INTRODUCTION 

Because of the difficulty of obtaining fully developed turbulent boundary -layer flows 
on flat plates and hollow cylinders at hypersonic speeds (ref. 1), turbulent boundary layers 

*This report was one of the papers presented at the Symposium on Compressible 
Turbulent Boundary Layers, Langley Research Center, Dec. 10-11, 1968. 



on nozzle walls have been surveyed extensively. Also, more accurate measurements can 
be obtained on tunnel -wall boundary layers because of their comparatively larger thick- 
ness. There have been only a limited number of comparisons made between the mea- 
sured profiles and integral thicknesses with the results of theoretical methods wherein 
the nozzle pressure-gradient history is taken into account. (See, for example, refs. 2, 

3, and 4.) Therefore, there is some question as to the magnitude of the effects of the 
favorable pressure-gradient history on the profiles measured at the nozzle exit. That 
is, because of past history, the boundary layer at the nozzle exit (where the pressure 
gradient is usually small locally) may not have the same characteristics as a boundary 
layer which developed in a zero-pressure-gradient situation. 

On the basis of skin -friction and heat-transfer data obtained on both a flat plate and 
a tunnel wall, Wallace (ref. 5) concluded that the two flows tend to have a close corre- 
spondence as far as surface phenomena are concerned. However, discussion in refer- 
ences 6, 7, and 8 (based on limited profile measurements) indicates that the temperature 
and velocity profiles for the nozzle-wall boundary layer are appreciably different from 
those for the flat-plate case. 

The purpose of this paper is to present some recent nozzle-wall turbulent- 
boundary-layer measurements at Mach numbers of 6, 8, and 19, and to compare the 
results obtained with two theories which utilize the upstream pressure history in com- 
puting the flow. One of the theories is of the conventional integral type, and the assump- 
tions used are examined as to relative applicability in flat-plate and nozzle flows. The 
other approach is a nonsimilar numerical solution of the governing differential equations 
and hence should be able to give an indication of the effects of the upstream favorable 
pressure gradient on the profiles at the nozzle exit. Numerical procedures for solving 
the differential equations have previously been developed by A. M. O. Smith (ref. 9) and 
others. 


SYMBOLS 

A constant in eddy viscosity expression (eq. (4)) 

Cp specific heat 

Cf skin -friction coefficient 

H total enthalpy 

H* form factor, b*/0 


2 


1 



l 

mixing length (eq. (4)) 


M 

Mach number 


N 

index in velocity profile (eq. (6)) 


P 

pressure 


Npr 

Prandtl number 


N Pr,t 

turbulent Prandtl number 


R 

Reynolds number 


r 

radius 


T 

temperature 


u,v 

velocity components in x- and y -directions, respectively 


x ,y 

Cartesian coordinates along and normal to surface, respectively 


a 

exponent in total-temperature-velocity relationship (eq. (5)) 


P 

density 


e 

eddy viscosity 


K 

eddy conductivity 


P 

viscosity 


T 

shear stress 


6 

boundary -layer thickness taken where — = 0.995 

u e 


6* 

displacement thickness 


e 

momentum thickness 

3 
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Subscripts: 
t total 

w wall 

00 local free stream 

e local external to boundary layer 

A prime denotes a fluctuating quantity. 

DESCRIPTION OF THEORIES 
Finite -Difference Approach 

The finite -difference method solves the compressible turbulent-boundary -layer 
equations in terms of mean flow quantities (ref. 10). With the following assumptions for 
the Reynolds stress and heat transfer 


-p u'v' = 

-p vT = 

and the definition of turbulent Prandtl number 

c p e 

= N 

K 

the governing equations become 
For continuity: 

i^( purj ) + J^( pvrJ ) = 0 (1) 

where j = 1 for axisymmetric flows and j = 0 for two-dimensional flows 


= K 


au 

3y 

9T 

9y 
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For conservation of momentum: 


Bu Bu 

pu — — + pv — — = 
9x By 


-& + + 
dx By M/By 


( 2 ) 


For conservation of energy: 


BH 8H 8 

p U -ii + p V — i = — 

8x By By 



1 + S. 


N Pr \BH L M e M 1 ~ N Pr,t ^ 

* N P r,tj9y ' V Pr M Pr N Pr>t j 


u 


Bu 

9y 


(3) 


The assumption is made that the eddy viscosity can be represented by the simple Prandtl 
mixing-length concept across the entire boundary layer. A Van Driest modification 
(ref. 11) is used near the wall, gas properties being evaluated at wall conditions. The 
expression for e is thus 


where 


e 




(4) 


A = 


26 ^L[ P JL\ /2 
Pvt \ T w J 


The variation of 1 / 5 with y/S was based on results from reference 12 and was 
assumed to be invariant with Mach number. The numerical values used in the present 
work are given in table I. For the calculations shown herein, Np r ^ was taken as 0.9 
(ref. 13) and was assumed to be constant. 

The method of solution of equations (1), (2), and (3) is similar to that used in refer- 
ence 14 and is essentially a linearized form of the implicit finite-difference numerical 
method of reference 15. As in reference 14, a variable step size in the y-direction is 
used. 


Integral Approach 


In the integral method the momentum integral equation for axisymmetric flow 
(eq. (42), ch. 9 of ref. 16) is solved by a variable -step -size fifth-order Runge-Kutta 
numerical scheme. In order to solve this equation, H* and Cf variations must be 
specified. In the outer part of the boundary layer, a temperature-velocity relationship 
of the form 


e s 


TfT w 




(5) 
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is assumed where the value of a depends on the flow configuration for which calcula- 
tions are being made. Further discussion of the values assigned to a is given in a 
subsequent section. Near the wall the total-temperature — velocity relation is obtained 
from the Colburn form of the Reynolds analogy. The wall total-temperature — velocity 
relationship is joined to equation (5) by an intermediate linear relationship which matches 

the wall equation at — = 0.01 and equation (5) at — = 0.10. 

u e u e 

The velocity profiles were assumed to be 



where N is a specified function of R e q and T W ^T^. For calculations of the nozzle- 
wall boundary layers, this function was taken as the solid line through the data shown in 
figure 1. The data shown are taken from references 4, 18 to 22. The wall -temperature 
factor used in the correlation is based on flat-plate type of data at the same Mach number 
and R e> q values but at different values. Since the present objective is the 

calculation of air and nitrogen axisymmetric nozzle-wall boundary layers, the data used 
to define an N variation are limited to these conditions. For those cases where only 
pitot data were available (shown flagged), the velocity profiles were obtained by using a 
value of 2 for a in equation (5). (See next section.) The velocity and temperature dis- 
tributions just described (eqs. (5) and (6)) were used to obtain H*. 

The skin-friction-coefficient expression of Spalding and Chi (ref. 17) was used 
directly, without any attempt to account for the effects of pressure gradient. 

RESULTS AND DISCUSSION 

Dependence of Temperature-Velocity Relation on Flow History 

As mentioned in the Introduction, the relation between total temperature and 
velocity may be different for flat-plate and nozzle-wall boundary layers. (See refs. 6, 

7, and 8.) A possible cause of this difference is the favorable pressure -gradient history 
of the nozzle flows. In order to investigate this difference further, all available data 
where heat transfer was present and the pressure gradient small locally, are shown 
plotted in figures 2 and 3 in the usual "Crocco" variables ^0 and jj-j. In figure 2 are 

shown the data which have been obtained on the flat-plate type of configurations, that is, 
flat plates, hollow cylinders, and in flows where the upstream pressure gradient was 
small (refs. 24 to 32). Included are unpublished measurements obtained on a flat plate 
at Mach 6.5 by Hopkins and Keener of NASA Ames Research Center. The amount of 
scatter exhibited by the data is appreciable, and this scatter may be due to the difficulty 
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of obtaining accurate measurements in the relatively thin boundary layers usually pres- 
ent in this class of configuration. Also, "flat-plate -flow" experiments do not always 
achieve the uniform conditions implied in the name. Effects due to induced pressures, 
wall -temperature variations, secondary flows, and proximity of the measuring station to 
transition may contribute to the scatter in the available data. However, the mean of the 
data is perhaps adequately represented for engineering purposes by the linear relation 
shown. As is well known, this linear variation is a solution to the governing differential 
equations for the special case of zero pressure gradient and turbulent Prandtl number 
equal to 1. Since the turbulent Prandtl number is probably of the order of 1 (ref. 13, for 
example), the "agreement" between the flat-plate data and the linear relationship seems 
to be reasonable. 

Shown in figure 3 are the available data (refs. 8, 20 to 22, and 33 to 37) from the 
"nozzle wall" class of configuration, that is, data obtained on wind-tunnel walls or flat 
plates and cones where a nozzle type of upstream pressure history was imposed. Included 
in this figure are the results of the present Mach 8 and 19 measurements. It is apparent 
that, in contrast to the data in figure 2, the nozzle data are better represented by a quad- 
ratic type of relationship (a = 2). R. K. Matthews of Arnold Engineering Development 
Center indicated that unpublished nozzle wall data at Mach 6, 8, and 10 from Arnold 
Engineering Development Center are also in general agreement with this quadratic trend, 
as are the data from reference 23. 

An indication of possible effects of the upstream pressure gradient history that 
might account for the apparent difference between the nozzle -wall and flat-plate 
temperature -velocity relation can be obtained from the present finite-difference calcula- 
tion method. Calculations were made for the boundary -layer development down the con- 
tour of the nozzle from which the Mach 8 data were obtained. The calculation was ini- 
tiated at the nozzle throat and carried through to the nozzle exit. The variation down the 
nozzle of the parameters needed for the calculation is shown in figure 4. 

Computed total temperature and velocity profiles are shown plotted in figures 5 
and 6 for several positions along the nozzle. The input profiles for both u/u e and 9 
at the throat were assumed to be 1/7 -power-law profiles, that is, the input total- 
temperature — velocity relation was linear. These profiles are indicated in figures 5 
and 6. Downstream of the throat in the favorable pressure-gradient region at 
x » 5 inches (12.7 cm) (station 2) the velocity profile has bulged outward, but the total- 
temperature profile has actually become less full. As can be seen from figure 6, these 
results lead to a temperature -velocity variation similar to that measured further down- 
stream on the nozzle walls. (See fig. 3.) At the measuring station near the nozzle exit, 
where the pressure gradient is locally small, the velocity profile is still considerably 
fuller than the total-temperature profile. This condition perhaps indicates that the 
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boundary layer has not yet relaxed to an equilibrium state. Profiles computed down- 
stream indicate that the temperature -velocity relation tends to approach the linear vari- 
ation with increasing downstream distance. (Further calculations have indicated that the 
predicted temperature -velocity relation at the nozzle exit is not a strong function of the 
input profiles at the throat.) 

There is another mechanism present in at least some of the available nozzle-wall 
investigations which could produce a decrease in the total temperature in the wall bound- 
ary layer. In the Mach 6 measurements reported herein, it was found that there was a 
layer of air adjacent to the settling chamber liner which was much cooler than the rest 
of the flow. If this thermal layer expands down the nozzle without mixing, the nozzle- 
wall boundary layer would be growing in a cooler effective free stream until the cool- 
layer mass flow is absorbed or swallowed in the wall boundary. Estimates of this 
"swallowing" distance give values of up to half the nozzle length for some cases, and it 
is possible that the effect within the boundary layer of the cooled mass could persist 
downstream much farther. 

The previous discussion has indicated that turbulent boundary layers on nozzle 
walls at positions where measurements are usually made may not be in equilibrium in 
the sense that the effects of upstream pressure gradients and/or thermal layers can per- 
sist for large downstream distances. However, the nozzle data are useful for at least 
two purposes: (1) The data can be used to improve prediction methods for boundary- 

layer displacement effects in hypersonic nozzle design; and (2) these data can provide 
valuable test cases for the development of nonsimilar finite -difference methods. That 
is, the numerical solutions for nozzle flows can be used to evaluate the range of condi- 
tions for which models of the turbulent shear term (usually based on equilibrium data) 
are applicable. 

Application of Theories to Mach 6, 8, and 19 Nozzle-Wall Data 

Mach 8 measurements at nozzle exit .- Pitot-pressure and total -temperature mea- 
surements have been obtained by William V. Feller and Robert A. Jones of Langley 
Research Center in the wall boundary layer at the nozzle exit of the Langley Mach 8 
variable -density wind tunnel which is described briefly in reference 38. These data in 
the form of velocity and Mach number profiles and plots of 6 against u/u e are shown 
in figures 7, 8, and 9. Predictions from the integral and finite-difference methods are 
also indicated. (A tabulation of the faired experimental profiles is given in table n.) 

The measured and theoretical velocity distributions are plotted in figure 7 as a 
function of y, the physical distance from the wall. Both methods predict the velocity 
profile data with acceptable engineering accuracy. As can be seen in the table included 
in the figure, both theories give a fair prediction of the integral parameters 9 and 6*. 
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The corresponding experimental Mach number profile shown in figure 8 is pre- 
dicted accurately by the integral theory over most of the boundary layer. The finite- 
difference theory underpredicts the Mach number through most of the boundary layer and 
indicates that the theoretical total -temperature profile is too full since the velocity pre- 
diction was correct. The total-temperature — velocity plots for the Reynolds number 
range of the tests are shown in figure 9. All the data are more or less in agreement 
with the quadratic type of temperature -velocity relationship. Also, within the data 
scatter there seems to be little effect of Reynolds number. The finite -difference theory 
predicts too rapid a recovery toward the linear variation (or a final total -temperature 
profile that is too full). This result is reasonable since the eddy viscosity model used 
(eq. (4)) is based on equilibrium data and hence might be expected to drive the flow 
toward equilibrium too rapidly. (Further evaluation of the present eddy viscosity model 
is presented in ref. 39.) However, the additional mechanism referred to previously 
(enthalpy deficit in settling chamber flow) could cause reduced total temperatures within 
the boundary layer and hence might be partly responsible for the disagreement between 
theory and data. 

Mach 19 measurements at nozzle exit .- Pitot -pressure and total -temperature mea- 
surements at Mach 19 have been obtained by William D. Harvey and Frank L. Clark of 
Langley Research Center in the wall boundary layer at the nozzle exit of the hypersonic 
nitrogen facility at the Langley Research Center which is described in reference 40. 
Velocity and Mach number profiles and total temperature -velocity plots obtained from 
these data are shown in figures 10, 11, and 12 along with the predictions of the integral 
method. (A tabulation of the faired experimental profiles is given in table III.) 

In the outer section of the velocity profile shown in figure 10, the prediction of the 
integral theory agrees with the measurements. However, in the inner region the predic- 
tion deviates from the data considerably. This discrepancy is due to the fact that a 
power-law representation of the profile is not valid in the viscous sublayer which, for 
the present conditions, extends out to values of y that are from 20 to 40 percent of the 
boundary -layer thickness. The accuracy of the present integral method is therefore 
limited by the assumption of a power -law velocity profile. The prediction of the integral 
theory for the Mach number profile shown in figure 11 is not very satisfactory. The rea- 
son for disagreement in the inner region of the boundary layer is probably the overpre- 
diction of the velocity as seen in figure 10. The variation of total temperature with veloc- 
ity indicated by the data is shown in figure 12. The data are seen to be above the 
quadratic variation. 

Mach 6 measurements at and downstream of nozzle exit .- Pitot-pressure and total- 
temperature measurements have been obtained by William V. Feller and Robert A.Jones 
of Langley Research Center in the boundary layer of the Langley Mach 6 high Reynolds 
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number tunnel both at the nozzle exit and at several stations along a straight-pipe exten- 
sion downstream of the nozzle exit. This facility has recently become operational and 
has a contoured axisymmetric nozzle with a 6° maximum turning angle and a nominal 
test-section diameter of 12 inches (30.48 cm). Velocity and Mach number profiles at the 
nozzle exit are shown in figures 13 and 14, whereas the temperature -velocity data for the 
various stations along the straight extension are shown in figure 15. 

The velocity profile data shown in figure 13 are in good agreement with the integral- 
theory prediction, as are the integral parameters 9 and 5*. (See table in fig. 13.) The 
predicted Mach number profile shown in figure 14 is also in reasonable agreement with 
the data. 

The variations of total temperature with velocity are shown in figure 15 for various 
positions along the straight -pipe extension. The data at the first three measuring stations 
(94 in. (238.7 cm), 124 in. (315 cm), and 172 in. (436.7 cm) downstream of the throat) are 
seen to be nearly the same and are in general agreement with quadratic variation. How- 
ever, the data at the last measuring station, which is some 60 boundary -layer thicknesses 
(using an average 5 of 2 in. (5.08 cm)) downstream of the nozzle exit, show a tendency 
towards the linear variation. In an unpublished U.S. Naval Ordnance Laboratory investi- 
gation by Lee and Yanta at some 30 boundary -layer thicknesses downstream of their noz- 
zle exit (start of dp/dx ~ 0 flow along straight section), they were unable to detect any 
tendency toward a linear variation. Thus large distances may be necessary in nozzle 
flows before recovery toward a linear temperature-velocity variation occurs. 

CONCLUDING REMARKS 

Comparison of available flat-plate and nozzle-wall turbulent-boundary-layer data 
indicates that the variation of total temperature with velocity is generally different in the 
two types of flow. Data for the flat-plate flows show considerable scatter but the average 
exhibits a linear type of relationship as opposed to the more quadratic variation always 
found in the nozzle data. Results of nonsimilar finite -difference calculations for the noz- 
zle flows indicate that at least part of this difference is caused by the favorable pressure- 
gradient history. 

Profile measurements along a straight section downstream of the nozzle exit of a 
Mach 6 tunnel indicate that a distance on the order of 60 boundary -layer thicknesses may 
be necessary before the flow begins to revert (or relax) toward a linear temperature- 
velocity variation which is perhaps typical of flat-plate type of flows. 
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The simple integral method used herein includes correlations taken from nozzle 
wall data and gives results in good agreement with nozzle profile measurements at 
Mach 6 and 8 and in reasonable agreement at Mach 18. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 18, 1969. 
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Figure 3.- Total temperature-velocity data for "nozzle wall" flows. 
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Figure 5.- Velocity and total temperature profiles predicted by finite-difference theory for various positions along Mach 8 nozzle. 
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Turbulent boundary-layer velocity profile at exit of Mach 8 contoured axisymmetric tunnel. R » = 3.1 x 10 4 ; M ( 
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Figure 10.- Measured turbulent boundary-layer velocity profile at exit of Mach 19 contoured axisymmetric tunnel compared with prediction of 

integral theory using a = 2.0; R g Q = 5.17 x 10 3 ; j^~= 0.17. 
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